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We study a class of models for brittle fracture: elastic theory models which allow for cracks but not 
for plastic flow. We show that these models exhibit, at all finite temperatures, a transition to fracture 
under applied load similar to that at a first order liquid-gas transition. We study this transition 
at low temperature for small tension. We discuss the appropriate thermodynamic limit of these 
theories: a large class of boundary conditions is identifled for which the energy release for a crack 
becomes independent of the macroscopic shape of the material. Using the complex variable method 
in a two-dimensional elastic theory, we prove that the energy release in an isotropically stretched 
material due to the creation of an arbitrary curvy cut is the same to cubic order as the energy release 
for the straight cut with the same end points. We find the normal modes and the energy spectrum 
for crack shape fluctuations and for crack surface phonons, under a uniform isotropic tension. For 
small uniform isotropic tension in two dimensions we calculate the essential singularity associated 
with fracturing the material in a saddle point approximation including quadratic fluctuations. This 
singularity determines the lifetime of the material (half-life for fracture), and also determines the 
asymptotic divergence of the high-order corrections to the zero temperature elastic coefficients. We 
calculate the asymptotic ratio of the high-order elastic coefficients of the inverse bulk modulus and 
argue that the result is unchanged by nonlinearities — the ratio of the high-order nonlinear terms 
are determined solely by the linear theory. 
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I. INTRODUCTION 

Early in the theory of fracture, Griffith used Inglis' 
stress analysis of an elliptical flaw in a linear elas- 
tic material to predict the critical stress under which a 
crack irreversibly grows, causing the material to fracture. 
Conversely, for a stressed solid the Griffith criterion de- 
termines the crack nucleation barrier: if the material has 
micro-cracks due to disorder or (less commonly) thermal 
fluctuations, how long does a micro-crack have to be to 
cause failure under a given load? In a sense, a solid un- 
der stretching is similar to a supercooled gas: the point 
of zero external stress plays the role of the liquid-gas con- 
densation point. Fisher's |^ theory of the condensation 
point predicts that the free energy of the system develops 
an essential singularity at the transition point. In this 
paper we develop a framework for the fleld-theoretical 
calculations of the thermodynamics of linear elastic the- 
ory with cracks (voids) that naturally incorporates the 
quadratic fluctuations, and we calculate the analogue of 
Fisher's essential singularity. As it is well known, the 
imaginary part of this essential singularity can be used 
to give the lifetime to fracture: what is the rate per unit 
volume of a micro-crack fluctuations large enough to nu- 
cleate failure? 

There is much work on thermal fluctuations leading to 
failure at rather high tensions, near the threshold for in- 
stability (the spinodal point) there is also work on 
the role of disorder in nucleating cracks at low tensions 
1^. We are primary interested in the thermal statistical 
mechanics of cracks under small tension. We must ad- 



mit and emphasize that, practically speaking, there are 
no thermal crack fluctuations under small tension — our 
calculations are of no practical significance. Why are we 
studying thermal cracks in this formal limit? First, for 
sufficiently small tension, the bulk of the material (ex- 
cluding regions near the crack tips) obeys linear elastic 
theory, thus making analytical analysis of the fracture 
thermodynamics tractable. Second, the real part of our 
essential singularity implies that nonlinear elastic theory 
is not convergent! Just as in quantum electrodynamics 
Q and other field theories |^], for all finite temperatures, 
nonlinear elastic theory is an asymptotic expansion, with 
zero radius of convergence at zero pressure. We will cal- 
culate the high-order terms in the perturbation expansion 
governing the response of a system to inflnitesimal ten- 
sion. We flnd it intriguing that Hook's law is actually a 
first term in the asymptotic series. 

The paper is organized as follows. In the next section, 
following methods known in the crack community ||^ — 
pT| , we carefully examine the thermodynamic limit of an 
equilibrium linear elastic theory with voids. We consider 
a crack as a special case of a void. We specify the class of 
boundary conditions which insure that the energy release 
is independent of the shape of the material boundary 
at infinity and independent of the prescribed boundary 
conditions. This is extremely important for the inves- 
tigation of the singular structure of the free energy, for 
the latter can develop singularities only in the thermody- 
namic limit [ p2[ . Using the complex variable method in 
a two-dimensional elastic theory, we calculate the energy 
release of an arbitrary curvy crack to quadratic order 
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in kink angles in section III. In section IV we find the 
spectrum and the normal modes of the boundary fluc- 
tuations (surface phonons) of a straight cut under uni- 
form isotropic tension at infinity. Section V is devoted to 
the calculation of the imaginary part of the free energy. 
The calculation of the contribution of thermal fluctua- 
tions depends on the "molecular structure" of our mate- 
rial at short length scales — in field theory language, it 
is regularization dependent. We calculate the imaginary 
part of the free energy both for C-function and a particu- 
lar lattice regularization, and determine the temperature 
dependent renormalization of the surface tension. Earlier 
we showed that the thermal instability of an elastic 
material with respect to fracture results in non-analytical 
behavior of the elastic constants (e.g. the bulk modulus) 
at zero a ppl ied stress. In section VI we extend the cal- 
culation [|l3| of the high order expansion of the inverse 
bulk modulus by including quadratic fluctuations. We 
show there that the asymptotic ratio of the high order 
elastic coefhcients, written in terms of the renormalized 
surface tension, is independent of regularization (for the 
cases we have studied), and we argue also that they are 
independent of nonlinear effects near the crack tips. (The 
asymptotic nonlinear coefficients depend only on the lin- 
ear elastic moduli.) In section VII we perform the sim- 
plified calculation (without fluctuations) in several more 
general contexts: anisotropic strain (nonlinear Young's 
modulus), cluster nucleation and dislocation nucleation, 
and three-dimensional brittle fracture. We also discuss 
the effects of vapor pressure — nonperturbative effects 
when bits detach from the crack! Finally, we summarize 
our results in section VIII. 



II. THE THERMODYNAMIC LIMIT OF THE 
ENERGY RELEASE 

Elastic materials under a stretching load can relieve 
deformation energy through the formation of cracks and 
voids. The famous Griffith criteria for a crack prop- 
agation is based on the balance between the energy re- 
lease and the increase in the material surface energy due 
to extending the crack. For a finite size system the en- 
ergy release is a well defined quantity that depends on 
the shape of the material boundary. The situation be- 
comes more subtle in case of an infinite elastic media. In 
principle one can calculate the energy release analyzing 
stress fields near the crack tips and thus avoid the neces- 
sity of worrying about infinite-sized media. This method, 
developed by Irvin in the 1950s, is known as the stress 
intensity approach Despite its enormous practical 

importance in numerical calculations, it is usually of little 
help in analytical calculations. To apply the stress inten- 
sity factor approach, one has to be able to compute the 
stresses near the crack tips, which is possible only in sev- 
eral simple cases. (The extension of the Irvin's method 
though the use of the path independent J-integral , 



for example, is applicable only for cracks with flat sur- 
faces.) Alternatively, the energy release can be calculated 
considering the system as a whole. In this approach, to 
compute the energy release one has to evaluate the work 
done by external forces and the change in the energy 
of elastic deformation. The change in the energy of the 
elastic deformation involves the difference between two 
infinitely large quantities for an infinite material; the lat- 
ter thus requires some sort of infinite- volume limit. In 
this section we discuss the energy release due to the re- 
laxation of the boundaries of a finite number of voids 
(cracks) in the limit of an infinitely large stressed mate- 
rial. The methods of this section will be used again later 
in the paper. 

We focus on the energy release calculation for a void 
formed in an infinite two-dimensional elastic material. 
The result is then extended to the case of a finite num- 
ber of voids (remember that a crack can be considered as 
a degenerate void) and for the energy of void formation in 
a three dimensional elastic material. The energy release 
for a finite size system with a crack has been calculated 
by Bueckner |^ and later generalized by Rice |^ for void 
formation. Their analysis allow for a large class of bound- 
ary conditions, mixing regions of fixed displacements and 
fixed stresses at the perimeter. In what follows we will 
use a slight modification of Bueckner's argument. 
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FIG. 1. To calculate the energy release in the thermody- 
namic limit we specify displacements along Ou and stresses 

along Oct- 

We define the energy release due to the formation of 
a void in an infinite material as the energy in a previ- 
ously stretched material, released by cutting out a hole 
in it and letting the hole boundary relax the stress: the 
deformation energy of the discarded piece is excluded. 
Let's consider an infinite linear elastic material subject 
to stress fields ct!^, and cr^ at infinity. We want to 
calculate the energy release E'leioasc due to cutting out a 
hole with boundary F/j. We assume that the hole bound- 
ary is stress free and non self-intersecting as a result of the 
stress relaxation. Fb denotes a regularization boundary 



2 



with a characteristic size L. Before the void formation, 
displacements along Tt are C/f = {u\,v\), while those 
along the prospective void contour Th correspondingly 
U'l — {u^jVi). Let Oct and Ou be nonintersecting parts 
of Fb, such that Oa U Ou = Tf,. Along Oa we fix the 
stresses cr^, cr^ and as a function of position; along 

Ojj we fix the displacements to be tj^ (Figure |^). In anal- 
ogy to the usage in field theory, each choice of boundary 
conditions O^, cr°°, Ojj, U\ we call a regularization. For 
a fixed L we calculate the energy release -Breioaso- The 
thermodynamic limit of the energy release is then given 
by 



lim ELca.sc- 



(1) 



Physically the described regularization means that we 
compute first the energy release due to the formation of 
the void Th in a finite size material with the boundary 
Ff,, and then push the outer (regularization) boundary to 
infinity. We will show that defined as above, i?reiease is 
independent of a particular choice of and Ou, even if 
they are themselves functions of L. The boundary condi- 
tion with Oa = Fb is known as a "fixed tension boundary 
condition" , while that with Ou = Fb is referred to as a 
"fixed grip boundary condition" . 

Let's first consider a straight cut. As shown by explicit 
calculation ||l^ , the "fixed tension" and the "fixed grip" 
boundary conditions for the special case of a straight cut 
of length i opened by a uniform isotropic tension T, give 
the same energy release 



32/i 



4y 



(2) 



The material elastic constants /i and x can be expressed 
through its Young's modulus Y and Poisson ratio a as 
follows 



Let e'",]\ (T,-)'' be the strain and stress fields of the first 



state before the cut is made and 



,(2) 



erf?'' define the 



fields of the second state with the void; finally displace- 
ments along Oct for the second state are U2 — {'^2' ^2) 
and displacements along the void boundary F;i are given 
by U2 — {u2,V2)- As the void boundary relaxes, the 
external forces do work 



We 



(5) 



where F„ — {Fn^ , Fn^ ) is the traction along Oa defined 
through the asymptotic stress fields 



Fn, = 



(6) 



with (n^^ny) being the outwards normal to Oa- The 
work done by the internal forces is the change in the en- 
ergy of the elastic deformation of the first and the second 

states 
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a-^-ei^ dA-- 



(7) 



The integration in (^ is performed over the area A of 
the material excluding the void; the summation over re- 
peating indexes is assumed. As a consequence of Hook's 
law, 



^^dA 



so we can rewrite (0) as 



(1) 



J2) 
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(1) 



= (2) 



)dA. 



(8) 



(9) 
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3-0- 



(3) 



X = 
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(The given value for x corresponds to a plain stress in 
a three dimensional elastic theory; for a plain strain one 
should use x — 3 — 4ct.) Note that (||) coincides with Grif- 
fith's result. This result is not trivial! If one calculates 
the energy change in a region Ft embedded in an infinite 
medium (i.e., neither fixing the displacements nor allow- 
ing them to relax at fixed stress) , the energy release does 
depend on the shape of the boundary. The relaxations at 
the boundary scale hke 1/L: even as L ^ cx), integrated 
over the perimeter they must be included for a sensible 
thermodynamic limit. 

For a general void, the energy release is a sum of the 
work performed by the external {We) and the internal 
(Wi) forces as a result of the void formation 



^release — 



W,. 



(4) 



Let's introduce a plus state specified by the stresses 



fields e 



(2) . 

and a minus state defined by the strain 



=(1) 



J2) 



Then the work of the internal 



forces (^ is a mixed energy of the plus and the minus 
states. According to Betti's theorem the mixed en- 
ergy equals one half the work done by the stresses of one 
state over the displacements of the other, no matter from 
what state the stresses or displacements are taken. With 
the stresses from the plus state and the displacements 
from the minus state we obtain 

W,,^l f F+U,-dh + l<f F+U^din (10) 



{F^+F,,){Ul-U^2)dh 
{F,,+Q){U^-U^)d£h 



FniU^~U^)di, 



Fh{Ui~U2)dih 
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where a traction of the first state Fh is defined as in m) 
through the stresses of the first state along F/j. From 
(H) and din) we find the regularized energy release M . 



release 



(11) 



The thermodynamic limit (^ is then taken by simply 
replacing U2 with its value 



- lim 

L—>-oo 



(12) 



for the infinite media. One can check (using for ex- 
ample the complex variable method described below) 
that the difference between two can at most be of order 
0{1/L), thus according to ( |Tl| ) and (|l]) not contributing 
to i?rcicasc- So wc coucludc that 



_ 1 

dease 



(13) 



Thus the total energy release is half the work done at the 
cut boundary. The above result is explicitly independent 
of a particular regularization from the discussed class as 
well as of the shape of the regularization contour F^. ft 
has a straightforward extension for a finite number of 
voids in a media. For a media with n holes with contours 



{F/j^ }, (13) generahzes to 



E, 



elease 
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(14) 



The same class of regularizations gives well defined ther- 
modynamic limit in a three dimensional case as well. The 
arguments are practically the same resulting in the en- 
ergy release (|l^) with the only change in the use of sur- 
face integrals rather than the contour ones for the regu- 
larization and void boundaries. 

In order to use ( p^ ) one has to know displacements 
along the void contour for the second state. For all but 
several simple cases this is a very complicated problem 
to approach analytically: at best, one can hope to get an 
asymptotic behavior of the displacement and stress fields. 
So for practical purposes we have to find an asymptotic 
analog of the above expression that would allow the cal- 
culation of the energy release from the asymptotic behav- 
ior of the displacement fields. To do this, it is convenient 
to specify the stresses at infinity — the regularization 
with Ojj = and = Ft,. As we have already shown, 
this does not restrict the applicability of the result. To 
get the asymptotic expression for the energy release we 
first note that elastic states one and two are each in equi- 
librium, so from Clapeyron's theorem p6| 



J2) (2) , . _ 



FnUiMb 



FnU^dh. 



FhUi'dih 



(15) 



Which, following (|), (|) and (|) gives 



E, 
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FnUldlb 



FhU^dih. 



(16) 



Next, we break U2 into two pieces: displacements of the 
second elastic state for the infinite media along F^, U^, 
and the boundary relaxation displacements, C/^, 



= Uf + Ur- 



(17) 



The boundary relaxation displacements comes from relax- 
ing the the stresses along Ff, for the infinite material to 
comply with the thermodynamic limit prescription: the 
"fixed tension boundary condition" . For the energy re- 
lease calculation it is sufficient to relax the stresses to 
0{1/L^) — the stresses of order 0{1/L^) generate the 
relaxation displacements along Ff, of order 0(1/L^) and 
thus do not contribute to the energy release in the limit 
L 00. So, from (§), (|l|) and (0) 

^release = ^^^^ ^ ^"(^/ + ^r)dib -\i F,,UlMb 



(18) 



This is the desired expression. One can see that (|^ is 
nothing but the difference between the regularized elastic 
deformation energy of the infinite media with the void 
(the first integral) and the elastic energy stored in the 
same media before the formation of void, not accounting 
for the elastic energy of the void itself (the remaining 
two integrals). Although the asymptotic expression for 
the energy release may look more complicated than ( [T^ ) 
and actually requires solving of two elastic problems, its 
use is justified because it is difficult to find displacements 
of the second state along the void contour. 



III. ENERGY RELEASE OF "SLIGHTLY" 
CURVY CUTS 

Armed with the knowledge that the energy release is 
independent of the shape of the regularization boundary, 
we now turn to the calculation of the energy release due 
to the opening of a "slightly" curvy cut F/,, in a two di- 
mensional isotropic linear elastic infinite media subject 
to a uniform isotropic tension T at infinity. We find an 
amazingly simple answer: the energy release of a curvy 
cut coincides to cubic order with the energy release of 
the straight cut with the same end points. 

An elastic state is completely defined once displace- 
ments (u, v) are known everywhere. Rather than consid- 
ering these two functions, Muskhelishvili ||l^ introduces 
two complex functions 4'{z) and 'ip{z) that in equilibrium 
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should be the functions of only one complex variable z 
(i.e. do not depend z). Moreover in our case (a uniform 
isotropic tension at infinity) ^{z) decomposes as 



00 (z) 



(19) 



The functions (/'o(z) and ''liiz) are holomorphic in the 
complex z plane including infinity but excluding the cut 
contour. This description associates the components of 
stress [iJxxi o'yy, o'xy) and displacement {u, v) to (0, tp) by 
the following relations 



(7xX+(7yy =2{(j)'{z) + (l)'{z)) 
7yy - (Jxx + '2,i(Txy = '^{Z(t)"{z) + i)' {x)) 

2fi{u + iv) — x4'{z) — zcf)'{z) — ^{z) 



(20) 



(The detailed discussion of the change of "variables" 
(m, v) (0, ^/j) along with the derivation of (p^)-([20[) can 
be found in [[l7[ .) Using a circular regularization contour 
(Ff, is a circle of radius L) and an expression analogous to 
(p^), Sih and Liebowitz explicitly computed stresses 
along the outer boundary Ff, and found the energy release 



-^release = - — ( 1 + x)R-e[yi] 

4m 



(21) 



with Hi being the residue of ''^{z) at infinity (the 1/z 
coefficient in the expansion about z = oo). Of course, to 
determine yi we still need to solve the asymptotics of the 
elasticity problem. 

To illustrate the correspondence between the energy 
release of a curvy cut and the straight one with the same 
end points, let's consider a rare example where it is pos- 
sible to find an exact analytical solution. Suppose a ma- 
terial with a "smile" cut - an arc of a circle ABC - of total 
arc length i (Figure |^) is subject to a uniform isotropic 
stretching T at infinity. Expanding the exact answer in 
llTf about z = oo, we find 



"87 6*2 (3 -cose 



(22) 



which according to (Ell) gives the energy release Exbc 



2«2 



■(1 + x) 



Ssin^ 61/2 
612 (3 -cos 61) ■ 



(23) 



On the other hand, for a straight cut AC of length 
{£/d)shiO, the holomorphic function ipAc{z) has asymp- 
totic behavior [ITtI 



, , , Ti^sin^e ( 1 
^Ac(.)^-^^ + 0(^ 

resulting in the energy release i?AC 



7rT2£2 
i2[i 



(1 + x) 



sm 



For small Q we find from ( |23| ) and ( |25| ) the advertised 
result: the energy release of ABC coincides with that of 
AC to cubic order in but not to quartic order! 



^r2^2 



Eac = 



32^ 
32n 



9^ 

X)(l-y 



779^ 
720 



+ 0(r) (26) 



(1 



^ + O(06' 
45 ^ ' 




FIG. 2. The "smile" -like cut ABC and the straight cut AC 
result in the same energy release to cubic order in 6. 

We proceed now with the general proof. First, an ar- 
bitrary cut is approximated by a finite number of line 
segments, parameterized by the kink angles — the 
angles between consecutive kinks. The exact shape of 
the cut is then restored as the length of each link goes to 
zero (as their number goes to infinity). The energy re- 
lease is evaluated to cubic order in the kink angles, e.g for 
the n - kink regularization, the energy release En{{ai}) 
for a curvy cut with a fixed separation ip between the 
endpoints is approximated as 



= 1 3 = 1 



n 71 n 



+ E E E Ef^a^a.a,,, + 0(af ) (27) 

2—1 j — 1 m—1 

where is the energy release for a straight cut of length 



£p and the coefficients E, 



(1) c.(2) 



and E. 



(3) 



-^Jm depend only 

(24) positions of the kinks along the cut. We claim 
that all coefficients up to cubic order are zero, and thus 
the energy of a curvy cut and the straight one with the 
same endpoints can differ only at 0{af). 

That £"1^^ and E^j}^ (in fact, all terms odd in the kink 

(25) angles) are zero follows from a symmetry argument: cuts 
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(having the same number of segments with the corre- 
sponding segments being of the same length) with kink 
angles {a^} and {— a.;} respectively, are mirror images to 
each other with respect to the first link. The boundary 
condition for our problem (a uniform tension at infinity) 
is reflection invariant, so 



En{{a,}) ^ En{{-a,}) 



(28) 



which requires that all energy release terms odd in the 

(2) 

kink angles vanish. To calculate E^j for a given pair of 
indexes we can put all kink angles to zero except for 
and aj, reducing the n-kink problem to a two-kink one. 
From now on we will consider only the two-kink problem 
to quadratic order in the kink angles. 




FIG. 3. The two-kink cut ABCD can be considered as a 
deformation of a straight cut AD. 

We choose the coordinate system XY in the complex 
z plane in such a way that the ends of the two-kink cut 
are on the X axis, symmetric with respect to the Y axis 
(Figure ^). Assuming a uniform isotropic tension T at 
infinity we rewrite (pi|), explicitly indicating the depen- 
dence of the energy release on the kink angles 

i^2(ai,a2) = -^(l + x)Re[yi(ai,a2)] (29) 

where yi{ai,a2) is 1/z coefficient in the expansion of 
the function ip^z) at infinity. As discussed earlier in the 
section, ip{z) is a holomorphic function in the complex z 
plane including infinity (the extended complex plane) but 
excluding the two-kink cut. The other function 4>{z) that 
is necessary for the specification of the equilibrium elastic 
state satisfies (19) with (j)o{z) holomorphic in the same 
region as tp{z)- The analytical functions 4>{z) and ^'{z) 
must provide a stress free cut boundary, which following 
[p7t can be expressed as 



{z) + z(j)'{z)+ij{z) 



(30) 



where f = + iFy is the complex analog of the force 
acting on the portion of the cut boundary between points 
W and X. 

It is important to note that any two pairs of functions 
((/)J(z), 'ip^{z)) and (0§(z), -0^(2)) that are holomorphic in 
the extended z plane excluding the same curvy cut, and 
which provide the stress free cut boundaries to 0{a^) 
(^), can differ only by 0{a^) everywhere: 



5il){z) 



^i(z) 



O(a^) 
= 0(a3). 



(31) 



This follows explicitly from Cauchy's theorem, but also 
follows from the elastic theory. Each pair {(t>Q{z) -f 
Tz/2,il^^{z)) or {(l)l{z)+Tz/2,ip'^{z)) defines the equilib- 
rium elastic state with stresses of order 0{a^) along the 
cut boundary and uniform isotropic stretching T at in- 
finity. So, {54>{z),5ip{z)) corresponds to the equilibrium 
state with the specified stresses of order 0{a'^) along the 
cut boundary and zero tension at infinity. Thus ( |3l| ) 
follows because the response to this force within linear 
elastic theory must be linear. The above argument guar- 
antees that once we find (f>o{z) and tp{z) that satisfy the 
discussed constraints to 0{a^), we can use them to cal- 
culate the energy release of the curvy cut to quadratic 
order. 

Let functions (j)^{z) and ^"^{z) define the equilibrium 
elastic state of a material with a straight cut AD subject 
to a uniform tension T at infinity. (f>o{z) — (f>''{z) — Tz/2 
and ip'^iz) should then be holomorphic in the extended 
complex z plane excluding the straight cut and should 
provide stress free boundaries along AD. Muskhelishvili 
finds [ITtI 



^(z)_T 



dz 



(32) 



d^^iz) _ T 



z£r> 



(z^-£//i)^z^-i^^/4 



w 



dz 



(To obtain (^^(z) and il^'^{z) we integrate (|32|); the arbi- 
trariness in the integration constants refiect the ambi- 
guity in the displacements up to a rigid motion of the 
material as a whole.) Note that (I>q{z) and tp'^lz) can be 
"made" holomorphic everywhere in the complex z plane 
excluding the two-kink cut ABCD, and thus can serve 
as a good starting point for the construction of 0o(-2) 
and "0(2). The process of an analytical continuation is 
demonstrated by Figure |[ (j}Q{z) (or equivalently il^^{z)) 
is holomorphic in the z plane excluding the straight cut 
AD (4a). Removing the region ABCD A (4b), we make 
it holomorphic in the complex plane excluding ABCDA. 
Now we analytically continue 4>q {z) from the link AD (4c) 
into the removed region (the continuation is possible ex- 
plicitly using (^2|)). The obtained function becomes holo- 
morphic everywhere in the complex z plane excluding the 
two- kink cut ABCD (4d), moreover the original function 
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and the one obtained through the analytical continuation 
coincide outside ABCDA. 
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FIG. 4. Functions <j)o{z) a-nd i^''{z) holomorphic in the com- 
plex z plane excluding the straight cut AD (a) , can be "made" 
holomorphic in the complex z plane excluding the two-kink 
cut ABCD (d). 



The idea of construction the holomorphic functions 
4>o{z) and is simple: we start with the functions 

4>q{z) and ip^iz) and calculate to quadratic order the 
stresses along the two-kink cut boundary ABCD un- 
der the analytical continuation as described by Figure 
4. The stresses along the curvy cut boundary (Figure 
5) are then compensated up to quadratic order in the 
kink angles by introducing counter-forces along the orig- 
inal (straight) cut, leading to corrected functions 6(j)'^{z) 
and Sip^iz), where (j){z) = (ji^iz) + 5(j)''{z) + 0{a^) and 
ip{z) = ijj {z) + &ijf{z) + 0{a^). For the calculation of the 
energy release (g9|) we need the real part of the residue of 
'ip{z) at infinity: wc will show that the residue of Sip'^^z) 
at infinity is zero and thus the residues of '(/'(z) and ip''{z) 
at z = oo are the same — which means that the energy 
release for the curvy cut ABCD is the same as that of for 
the straight cut AD. 





FIG. 5. The stress free boundary of a two-kink ABCD cut 
(b) can be mimicked by applying the tangential force to the 
previously unstressed (a) straight cut boundary AD. 

Let's assume that points W and X are on the upper 
boundary of the link AB. From Figure ^, z = i -I- ip(t + 
ip/2) + 0{a^), where t e AB' and /3 = (1 - fci)ai + (1 - 
k2)a2 + 0(a3). Using iM) we find 



X 



w 



X 



11+ 



+ 0(a=^) 



w 



20\t-rlpl2f<\>\tr 



+ 0(a3), 



(33) 



w 



where t runs along AB'; the -I- superscript means that 
the values of (\>%(t) and '>\}[t) should be taken at the up- 
per boundary of the straight cut. To obtain the second 
expression in ( |33| ) one can plug in the explicit form (|3^), 
or - more elegantly - note that for t G AB', 0^(t) is pure 
imaginary and ip''{zy = —z(t)^{zy' . Either way, it follows 
that the functions 5(j)'^{z) and 5^'^{z) satisfy 



iSf 



S(f>''{t)+ + t6<f)^{ty' +6ip'={t) 



/+ 



w 



-2f3''{t + £p/2f^%ty'^ 



X 



w 



(34) 



This is the force we need to add along the straight cut 
just below segment AB to cancel the stress along the 
curvy cut. Similar expressions can be found for the forces 
needed below BC and CD. To find (50^(z) and 6i/j'^{z) we 
have to solve the elasticity problem for the material with 
the straight cut AD, subject to these applied forces iSf 
along the cut boundary. Fortunately, this problem allows 
a closed analytical solution m^. Expanding the exact 



expression for 5ip'^{z) in |17| we find 



Sil^^iz) 



Airiz 



Re 



iSfixia)) 



da + O 



(35) 



where the integration is along the unit circle 7 in the 
complex plane, and iSf is a function of a variable point 
x{a) = £p{a + along the straight cut boundary 

AD. Notice from (|34| ) that iSf is pure imaginary evalu- 
ated on the upper boundary of the link AB: in this case 
\t\ < £p/2, and so the argument of the square root in ( |3^ ) 
is negative resulting in pure imaginary ip'^ (t) , thus (t) 
is also pure imaginary. In fact, as it can be checked ex- 



plicitly. Re 



for arbitrary W and X along the 



cut boundary. So we conclude from ( pSD that the residue 
of Sip'^lz) at infinity is zero and thus the energy release 
for the curvy cut ABCD is the same as for the straight 
cut AD. The underling physical reason for this seeming 
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remarkable coincidence is that to imitate the stress free 
curvy cut to quadratic order in the kink angles we have to 
apply only tangential force along the straight cut (pure 
imaginary iSf means Fy = 0), which do no work because 
a straight cut under a uniform isotropic tension at infin- 
ity opens up but does not shrink 

We find that the energy release E{£p) of the curvy cut 
with projected distance £p between the endpoints is the 
same to quadratic order in the kink angles as the energy 
release of the straight cut of length £p. The latter one is 



given by the second formula in (2( 
coincides with Griffith's result (E)] 



j) with = (it also 



Eiip) 



32^ ^ 



(36) 



The natural variables to describe the curvy cut are its 
total length £ and its curvature k{x), x G [0,£]. In what 
follows we express (|3^) in these variables and find the 
normal modes of the curvature that diagonalize the en- 
ergy release. 

For the two-kink cut ABCD (Figure ^ of total length 
£ one can find 



/ — 



xi{£-xi) 2 
^'^ 2^"i 



xi{£ - X2) 



£2 



-aia2 



X2(£ - X2) 
2^2 ' 



(37) 



where xi and X2 parameterize the kink positions: the 
length of the link AB is assumed to be xi and the length 
of the segment ABC equals X2- Similarly, for the n-kink 
cut of total length £ with the kink angles {a^} parame- 
terized by their distance {xi} from the cut end 



1 



XiXj imn{xi,Xj) 



(38) 



Expressing the kink angles through the local curvature 
of the curve, at = k{xi)Axi/X, we find the continuous 
limit of (M) 



Jo 



^k{x)Mix,y)kiy) 



with 



~0{k{xf), 



M{x,v) 



xy min(a;, y) 
"P^ £ 



(39) 



(40) 



and the scale A is introduced to make the curvature di- 
mensionless (A can be associated with the ultraviolet cut- 
off of the theory — roughly the interatomic distance). 
Substituting (^0|) into ( |3q ) we find the energy release 
E{£,k{x)) — E{£p) of the curvy cut in its intrinsic vari- 
ables 



E{£,k{x)) = 
7rT2^2 
32^ 



(1 + x) 1 



Jo 



^k{x)M{x,y)k{y) 



+o{k{xY 



(41) 



To find the normal modes of the curvature we have to find 
the eigenvalues and eigenvectors of the operator M{x, y). 
If kn{x) is an eigenvector of M{x,y) with eigenvalue A„, 
then 



^M{x,y)k4y). 



(42) 



From (B^), M{0,y) = and M{£,y) = for arbitrary 
y € [0,F[7so from ( ^ the eigenvectors of M{x,y) must 
be zero at a; = and x = £: fc„(0) = kn{£) = 0. An 
arbitrary function kn{x) with this property is given by 
the Fourier series 



m—1 



2A irmx 
£ ^^"^ 



(43) 



where the overall constant ^2X/£ is introduced to nor- 
malize the Fourier modes with the integration measure 
dx/X over x S [0,^]. One can explicitly check from 
( ^ ) that each Fourier mode ^/2XJ£ sin{'!Tmx / £) is in 
fact an eigenvector of M{x, y) with the eigenvalue Am = 
£/ (TT'^m?X). In terms of the amplitudes of the normal 
modes {c„}, (|4l]) is rewritten as 



E{£ACn}) 



7:T' 



'21)2 



32/u 



-(1+x) 1-E 



+ 0{cl). 
(44) 



( ^ ) is the main results of the section: we've calculated 
the energy release of an arbitrary curvy cut in its in- 
trinsic variables — the total length £ and the curvature 
k{x), X G [Q,£] — to quadratic order in k{x) and found 
the normal modes of the curvature that diagonalize the 
energy release. 

In conclusion we mention that the measure in the kink 
angle space is Cartesian — fl - dui — (and thus the func- 
tional measure Dk{x) = Da{x) is Cartesian), so the mea- 
sure in the vector space of the amplitudes of the normal 
modes {c„} is also Cartesian — Yi'^^i^'^n, because the 
Fourier transformation {fc(x)} —^ {c„} is orthonormal. 
This will be important in section V, where we will be 
integrating over crack shapes. 



IV. SURFACE PHONONS 

In the previous sections we have extensively discussed 
the calculation of the energy release due to the equi- 
librium opening of a cut in an elastic material. Since 
our goal is to deal with cracks as thermal fluctuations. 
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we must also deal with the more traditional elastic fluc- 
tuations — phonons, or sound. We find here that the 
bulk fiuctuations decouple from the new surface phonon 
modes introduced by the cut. We discuss the quadratic 
fluctuations for linear elastic material with a straight cut 
of length i subject to a uniform isotropic tension T at 
infinity; more specifically, we calculate the energy release 
for the material with an arbitrary opening of the straight 
cut and we find collective coordinates (normal modes) 
that diagonalize the change in the energy. 

An elastic state of the material can be defined through 
the specification of its displacements U = {u, v) at every 
point {x, y). For the material with a cut, the fields u{x, y) 
and v{x,y) can in principle have a discontinuity along 
the cut: assuming that the cut is an interval (x, y) = 
(H/2,£/2],0), 



2gx{x) — u{x,0+) — u{x,0~), for x £ [- 



-£/2,£/2] 
(45) 



and 



2gy{x) =v{x,0+) ~v{x,0~), for xe[~£/2,e/2] 



(46) 



may be nonzero. It is clear that the arbitrary state U 
can be decomposed into the superposition of two states 
Ug = {ug,Vg) ond Uc = (uc,Wc), whcrc Ug is the equilib- 
rium state for given displacement discontinuity {gx,gy) 
at the cut boundary (]45[-^6|) and tension T at infin- 
ity that maximizes the energy release, and Uc given by 
(uc,Vc) — (u — Ug,v — Vg) is a continuous displacement 
field everywhere. Recall that the energy release is the 
sum of the work done by the external forces and the 
work done by the internal forces (^). We define the en- 
ergy release E for the elastic state U with respect to the 
equilibrium state of the material Uq = (uo,vq) without 
the cut under the same loading at infinity, as a limit of 
this difference for finite size samples with boundary Tb 
and enclosed area A. We find following (|), (||) and (|) 



E= (t Tn{U - UQ)dl + - 



)dA 
(47) 



where ct,° and e?, are the stresses and strains of the equi- 

librium elastic state of the uncracked material Uq] aij 
and Cij are the stresses and strains of the elastic state 
of material with the straight cut and displacement field 
U; and H is a unit normal pointing outwards from the 
regularization boundary Th. This argument is similar to 
that in the second section, but the elastic state U is not 
an equilibrium one and so the arguments there are not 
directly applicable. We rewrite the energy release ( p^ ) 
making use of the decomposition U = Ug + Uc io get 



E = 



Tn{Ug - Ua)de + - 



<j!^e%)dA 



cr^jB^jdA 



TnUrde - l 



e7.+^^A^dA. 



1-3 V 



(48) 



The first two integrals in (|4^) give the energy release 
for the equilibrium elastic state with the specified cut 
opening {gx,gy) and tension T at infinity. According to 
our decomposition this energy release is maximum for 
given gx{x) and gy{x) and thus can not increase linearly 
by tuning Uc- The latter is true only if the last two 
integrals on the RHS of (^8|) - linear in Uc - cancel each 
other. (This can be verified explicitly by integrating by 
parts the last integral on the RHS of (|4^) and using the 
fact that Ug is an equilibrium state.) Thus 



E= i, Tn{Ug - Uo)de + - 



2 11. '^^J^v' 



(49) 



the energy factors, and the last term representing the 
continuous degrees of freedom does not "feel" the pres- 
ence of the cut and thus will have exactly the same spec- 
trum as that of the uncracked material. 

Although the elastic state Ug is an equilibrium one, the 
cut boundary is in general stressed, and so we still have 
to modify the result of the second section for the energy 
release ([T^). From the first equation in (p^), the elastic 
energy of the uncracked material is given by 



TnUode. 



(50) 



The elastic energy of the material with the cut (the 
second equation in (^5|)) is modified to incorporate the 
stressed cut boundary 



TnUgde - 



FhUgdl. (51) 



The second integral in (|51| ) is over the cut boundary Th 
and Fh is the force we have to apply to the cut boundary 
to insure its displacements satisfy (|4^p6|). With this 
change, following (|4g|-pT|) we find that the energy release 
as given by (Eq) decreases by 5E 



SE = 



1 



FhU.dL 



(52) 



In the spirit of the third section, the equilibrium elastic 
state Ug can be described by the analytical functions (j){z) 
and ^{z); (j^oiz) = (j){z) — Tz/2 and ^{z) are holomorphic 
in the extended complex z plane excluding the straight 
cut and are constrained to provide displacement discon- 
tinuity {gx,gy)- The energy release Eg is then smaller 
than the one given by (0^) by 6E 
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E, = -^{i + x)Myl]~5E 



(53) 



where is the 1/z coefBcient m the expansion of ^{z) 
at infinity. 

To determine the functions (t>a{z) and 4'{z) we con- 
formally map the complex z plane with the cut to the 
outside of the unit circle 7 (Figure ||) 



(54) 



so that the unit circle in the C plane is mapped to the 
straight cut boundary Th in the original plane, Th = 
w(7). 




FIG. 6. The determination of the holomorphic functions 
describing the equilibrium elastic state of the material with a 
straight cut is simplified in the conformal plane ^, where the 
unit circle 7 corresponds to cut boundary Th in the original 
z plane. 

The elasticity problem is reformulated in the confor- 
mal plane as follows: we have to find analytical func- 
tions 4>g{() = 4>{iv{()) and ipgiC) = ^{uj{C)), such that 
'^so(0 = 4'giC) ~ and V's(C) are holomorphic in the 
extended complex C, plane outside the unit circle and give 
the maximum energy release with displacement disconti- 
nuity {gx,gy) at the cut boundary. We introduce 



g{a) = Ug + ivg 



(55) 



where a = exp(iQ;), a e [0, 27r), is a parameterization 
of the unit circle 7. Since a and l/a represent opposite 
points across the cut, (BH-Esh require 



g{a)-g{l/<7) = 2[g4uj{a)) + igy{uj{a))], a € [0,^). 

(56) 

It is important to note that the equilibrium elastic state 
that maximizes the energy release for given displacement 
discontinuity (gx^gy) is unique; on the other hand ( |5^ ) 
determines only the asymmetric modes g'^y'^{<T) of the 
crack opening displacement for this state 

25-y-(a) = [ug{a) + iVg{a)] - K(l/a) + tVgil/a)] (57) 
= 2[gx{tj{a)) + igy{Lj{a))]. 

The symmetric modes g^^™{(j) 

2g^y^{<j) = [ug{a) + ivg{a)] + [ug{l/<j) + tVg{l/a)] (58) 



left unconstrained by (^6|), should then be relaxed to 
provide the maximum energy release for given g^^'^"^{a). 
Thus, to calculate the energy release Eg of the elas- 
tic state Ug we first find the energy release E{g) — 
^^(^asym _j_ ^sym-j ^^iQ equilibrium state with an ar- 
bitrary displacement along the cut boundary g{a) and 
then maximize the result with respect to g^^"^ 



E„ 



(59) 



In what follows we will use (fic.iC) and i'c^iC) to describe 
the equilibrium elastic state with an arbitrary displace- 
ment g{a) along the cut boundary and tension T at in- 
finity. (The energy release for arbitrary g{a) is still given 
by ( |52| - ^3|) .) Making the change of variables z w{C,) in 
( pO| ) and putting ^ = ct we obtain a constraint on (f>(;{C) 
and V'c(C) that guarantees the displacements along 7 to 
be g{a) 



Lo'{a) 



^c(c^) = 2^.9(0-) ■ 



(60) 



Once the solution (0^,1/),^) of the elasticity problem is 
found, we can compute the correction (^2|) to the energy 
release. Introducing the polar coordinates (p, 9) (Figure 



in the complex C, plane, Fh = Fpp + FgO and Ug 



Vpp + vg9, and using d£ — {10' {a)\\da\ we find 



(61) 



1 



o-ppWp + ap0V0 I |w'((T)||dcr| 



where in the second equality we express the force through 
the stress tensor components: Fp — Upp and Fg — Opo. 
The stress tensor components Opp and OpQ are given in 
terms of the 4>q{C) and ipciC) functions that, as we already 
mentioned, completely determine the equilibrium elastic 
state. Muskhelishvili finds IT^ 



app - lape - + == - ^==<i -fTT^MO 



(62) 



Noting that the transformation of the displacements 
along the unit circle in the Cartesian coordinates {ug, Vg), 
g ~ Ug + ivg, to the polar coordinates {vp, vg) is [|l^ 



1 uj'{a) 



Vp + ivg^ - \)" {A ug + ivg) ^ - ^ ff(cr), (63) 



we conclude from (O) 



1, 



6E = -Rej){app - iapg){vp + ivg)\uj' {<7)\\d(T\ (64) 



-Re * (cTpp - 1(7 pg) g{cr) — 

2 7-y a la 
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is given by (p2h with ^ a {p ^ 1) 



where a^n — i 

From ( |53f ) and (64) we find the energy release E{g) 



(i + x)Re[yi(g)] 



ttT 
4^ 



-Re (f) [(jpp - tape) .9(0-) — 



(65) 



where yi{g) is the 1/z coefficient in the expansion of 
'ip(;{oj~^ (z)) at z = CXI. 

The equilibrium elastic problem for material with the 
straight cut allows a closed analytical solution for the 
arbitrary specified displacement g{a) along the unit circle 
in the conformal plane Using the fact that 0Cq(C) and 
ip(^{C) are holomorphic functions outside the unit circle 
that satisfy (pOh, Muskhehshvih finds Il7| 



MO = 
MO = 



TiC 2fi 1 / 9icr)da T£ 



X '2.711 (7 - C 8xC 
W)d<J T£fx 2C 
a-C ^ 8 U C'-l 



(66) 



Assuming that g{a) is smooth, we represent it by a con- 
vergent Fourier series 



+00 



5(0') = ^{(^71 + ibn)<j''' 



(67) 



Using representation (|67| ) for g{a) we find from ( pq ) the 
energy release E 



E{g) = \ ^ a_i+xai 

+ 00 



(68) 



7rr^£^(l + x) / . , 1 
128^1 ^x 



(The computations are tedious, but straightforward: first 
we substitute ( |67| ) into ( |6^ ) to find the solution of the 
elasticity problem in terms of the Fourier amplitudes 
{on, 6„}, then we calculate the stress tensor components 
at the unit circle using (|62|) , and finally plugging the re- 
sult into ( p5| ) we obtain ^pq).) The next step is to relax 
the symmetric modes in the crack opening displacement 
given by 5(17). From ( |67| ) and ( ]5^ ) we find 



Ug — (a„ + a_„) cosna + (6_n — fe„) sinna (69) 

n=l 

+00 

) COS na + {an — a-n) sin na 



which with the change of variables 



(70) 



is rewritten as 



+00 

E 

n=l 

+00 



Vn cos na -|- Un sin na 



(71) 



iig — u„ COS na -I- u„ sin na. 



It is clear now that the asymmetric modes of the crack 
opening displacement are described by {u„,i;„}, while 
the symmetric ones are specified by {un,Vn\- (Recall 
that points parameterized by a and 1 /a (or equivalently 
a and —a) are opposite from one another across the cut.) 
The amplitudes {un, Vri \ are uniquely determined for the 
given [gx ,gy)- From (^^ and ( |7l| ) 

gx{^/2, cosa) -I- igy{t/2 cosa) = 2.{un + ivn) sin na 



(72) 

where a € [0, tt]. Using the transformation inverse to 
( [70I ) we can express the energy release (68) in terms of 
{w„, w„, u„, The obtained expression is maximum 

for 



X 
X 



(73) 



1 



Vi=Vi- 



1 + X 

1-X , Tt{x-l) 



1 + X 

and gives the energy release Eg 

+00 



8fi 



En 



Tin 



-Vl 



2TTfi / 2 I 2 



1 ^ 

n=l 



2 X 

Finally, the maximum of ([74|) is achieved for 

v^'''' = 

<- = 0, n^l 

„max _ 7rr^(l + X) 
- 8^^ 



(74) 



(75) 



and 



771 max 

^9 



nTH\l + x) 
32/x 



(76) 



n=l 
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which, as one might expect, corresponds to the equiHb- 
rium opening of the cut |17| and the energy release as- 
sociated with this opening (^. Expanding (74) about 



find 



we 



32/^ 



1 + X 



+ 00 

n=l 



n u„ + vt 



(77) 



Expression ( |77| ) is the desired result: we find that the 
crack opening displacements (specified on the unit circle 
in the conformal plane) 



r 1 J 1 ^ ^ 

\u, v\ — < Vnz — : — COS na + it„ smna. 



1 + X 

'■i + x 



COS na + Vn sm na 



imposed on the saddle point cut opening 



} = {0, 



8n 



4 



(78) 



(79) 



diagonalize the energy release and thus are the normal 
modes; with the excitation of the n-th normal mode with 
the amplitude {u„,ti„} the energy release decreases by 

2TTfin{ul + vp/{l + x)- 

Although ( |74|) has been derived for the material under 
uniform isotropic stretching at infinity, it can be rein- 
terpreted to describe the minimum increase in the en- 
ergy AE of the material under a uniform isotropic com- 
pression (pressure) P at infinity, due to the opening of 
the straight cut with specified displacement discontinuity 
along its boundary. For the dis pla cement discontinuity 
given by ( [72|) we find similar to (fz^ ) 



AE = vi 



27r^ 
X 



+ 00 

n=l 



n w„ -f vt 



(80) 



One can use the same arguments that lead to (^9|) to 
show that the crack opening normal modes ( |78| ) decou- 
ple from all continuous modes (that are present in the 
uncracked material) and thus leave their spectrum un- 
changed. The saddle point is however unphysical in this 
case: as follows from (|7^), (T in ( [79| ) should be replaced 
with —P), it corresponds to a configuration where the 
material overlaps itself. 



V. THE IMAGINARY PART OF THE 
PARTITION FUNCTION 

Elastic materials at finite temperature undergo a phase 
transition to fracture at zero applied stress, similar to the 
first order phase transition in spin systems below the crit- 
ical temperature at zero magnetic field. The free energy 
of an elastic material under a stretching load develops 



an imaginary part which determines the material life- 
time with respect to fracture. The imaginary part of the 
free energy has an essential singularity at zero applied 
stress. In this section we calculate this singularity at low 
temperatures in a saddle point approximation including 
quadratic fiuctuations. 

Consider an infinite two-dimensional elastic material 
subject to a uniform isotropic stretching T at infinity. 
Creation of a straight cut of length £ will increase the 
energy by 2a£, where a is the surface tension (the energy 
per unit length of edge), with a factor of 2 because of the 
two free surfaces. On the other hand, the cut will open up 
because of elastic relaxation. Using (|7^) for the energy 
release we find the total energy E{i) of the straight cut 
in equilibrium under stretching tension T: 



E{i) = 2ae - 



Introducing 



4 = 



32^ 



32/iQ; 



nT^il + x) 



we can rewrite the energy of the crack as 



E{e) = 2ae-a — 



(81) 



(82) 



(83) 



It follows that cracks with £ > £c will grow, giving rise to 
the fracture of the material, while those with £ < £c will 
heal — a result first obtained by Griffith |^]. At finite 
temperature a crack of any size can appear as a thermal 
fluctuation, which means that for arbitrary small stretch- 
ing T the true ground state of the system is fractured into 
pieces and so the free energy of the material cannot be 
analytical at T = 0. Because the energy E{£c) — a£c 
grows as 1/T^ as r — > 0, interactions between thermally 
nucleated cracks are unimportant at small T and low 
temperatures (allowing us to use the "dilute gas approx- 
imation"). 

The thermodynamic properties of a macroscopic sys- 
tem can be obtained from its partition function Z: 



Z=J2Y. exp(-/3S„Ar) 



(84) 



W=0 n 



where the summation N is over all possible numbers of 
particles (cracks in our case) and the summation n is over 
all states of the system with N cracks. 

To begin with, let's consider the partition function of 
the material with one cut Zi 



Zi = ^exp(-/?i?) 



(85) 



where the summation is over all energy states of the ma- 
terial with a single cut. The calculation of the imaginary 
part of the partition function is dominated by a saddle 
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point, that in our case is a straight cut of length l^. The 
straight cut is the saddle point because it gains the most 
elastic relaxation energy for a given number of broken 
bonds (we explicitly show in section III that curving a 
cut reduces the energy release). For now we neglect all 
fluctuations of the critical droplet (the cut of length £c) 
except for its uniform contraction or expansion — fluctu- 
ations in the length of the straight cut. Introducing the 
deviation A£ in the cut length from the critical length i^, 
M = e-L,vfe find from M^) 



(86) 



The fact that this degree of freedom has negative eigen- 
value means that direct computation of the partition 
function yields a divergent result. A similar problem for 
the three-dimensional Ising model was solved by Langer 
[^8| : one has to compute the partition function in a sta- 
ble state P = —T (compression), and then do an ana- 
lytical continuation in parameter space to the state of 
interest. The free energy develops an imaginary part in 
the unstable state, related to the decay rate for fracture 
p9[ : the situation is similar to that of barrier tunneling 
in quantum mechanics pO| , where the imaginary part in 
the energy gives the decay rate of a resonance. We have 
explicitly implemented this prescription for the simplified 
calculation of the imaginary part of the free energy : 
for the elastic material under a uniform isotropic com- 
pression at infinity allowing for the nucleation of straight 
cuts of an arbitrary length with an arbitrary elliptical 
opening (mode vi in (^0|)), we calculated the free energy 
in a dilute gas approximation. We carefully performed 
the analytical continuation to the metastable state de- 
scribing the elastic material under the uniform isotropic 
stretching T at infinity and found the imaginary part of 
the free energy 



A 



exp 



^r2(x + i) 



(87) 



where A is the area of the material and A is the ultravio- 
let cutoff of the theory. (The version of equation ( ^ ) as 
derived in |p^, overcounts the contribution from zero- 
restoring-force modes {27rA/X^) by factor 2. Because 
cracks tilted by 9 and n + 9 are identical, the proper 
contribution from rotations must be tt, rather than 27r.) 

The alternative to this analytical continuation ap- 
proach is to deform the integration contour over the am- 
plitude of the unstable (negative eigenvalue) mode from 
the saddle point A£ — along the path of the steep- 
est descent Q . More precisely, we regularize the direct 
expression for the partition function 



Zi = ZqItt 



A 



dM 



exp< — /? a£, 



a- 



(88) 



(which diverges at big ) by bending the integration 
contour from the saddle into the complex plane: 



Zi = Zo vr 



A\ f° dM 



A2 



^ exp<^ -/3( q;4 - a ^ 



(89) 



('^^) °^P(~/^"^c) 



dM r„ A^2 

_ exp /Ja- 



in (Hfg]) the factor (ttA/A^) comes from the zero- 
restoring-force modes for rotating and translating the 
cut, and Zq is the partition function for the uncracked 
material (unity for the present simplified calculation). 
The second integral in (|89| ) generates the imaginary part 
of the partition function 



\ f A\ f 7:1 

\.^Z, = ±-Zo (^TT- j exp(-/3a4) 



1/2 



(90) 



with the ± sign corresponding to the analytical contin- 
uation to either side of the branch cut of the partition 
function. (We showed in ||l^ that partition function is 
an analytical function in complex T with a branch cut 
along the line T e [0, +oo).) In a dilute gas approxima- 
tion the partition function for the material with N cuts 
Zjv is given by 



Z 



N 



Zo 



jZi/Zp) 
Nl 



N 



(91) 



which from (p4|) determines the material free energy 



F 



■InZ 



N 



InZo 



1^ 
PZo' 



(92) 



Following ( pO| ) and (£^) we find the imaginary part of the 
free energy 



Imf^^'-P'<=(T) = ±-|-('^ 



2\ 1/2 



/32TA2 V X + 1 



A 



(93) 



( |93| ) differs from (87) only because for the calculation 
of the imaginary part of the free energy in fl^ we used 
two degrees of freedom: the length of the cut and its 
elliptical opening, while in current calculation there is 
only one degree of freedom. One can immediately re- 
store ( |87| ) by adding the vi mode of (77) to the en- 
ergy of the elastic material ( p6| ) and integrating it out. 
From (fz?!), the vi mode generates an additional multi- 
plicative contribution Z^-^ to the partition function for a 
single crack Zi, and thus from (^2|) changes the imagi- 
nary part of the free energy for multiple cracks i^'simpic^ 

IniJ?simplc 



+ 00 



dvi 



■ exp 



2nfil3 



1/2 



2/i/3A2 



(94) 
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which will cure the discrepancy between ( P3| ) and (87). 
Although the analytical continuation method is theoret- 
ically more appealing, the calculation of the imaginary 
part through the deformation of the integration contour 
of the unstable mode is more convenient once we include 
the quadratic fluctuations. It is clear that both methods 
(properly implemented) must give the same results. 

We have already emphasized that the above calcula- 
tion ignores the quadratic fluctuations about the saddle 
point (except for the uniform contraction or extension of 
the critical droplet), which may change the prefactor in 
the expression (|9^) for the imaginary part of the free en- 
ergy and may renormalize the surface tension a. There 
are three kinds of quadratic fluctuations we have to deal 
with. (I) Curvy cuts — changes in the shape of the tear 
in the material: deviations of the broken bonds from a 
straight-line configuration. (II) Surface phonons — ther- 
mal fluctuations of the free surface of the crack about its 
equilibrium opening. (Ill) Bulk phonons — thermal fluc- 
tuations of the elastic media that are continuous at the 
cut boundary. To incorporate these fluctuations we have 
to integrate out the quadratic deviation from the sad- 
dle point energy coming from their degrees of freedom 
(as we did for the surface phonon vi above). In all cases 
the answer will depend upon the microscopic lattice-scale 
structure of the material. In field-theory language, our 
theory needs regularization: we must decide exactly how 
to introduce the ultraviolet cut-off A. Here we discuss 
the lattice regularization, where the cut-off is explicitly 
introduced by the interatomic distance, and ^-function 
regularization, common in field theory. We find that the 
precise form of the surface tension renormalization and 
the prefactor in the imaginary part of the free energy 
depends on the regularization prescription, but certain 
important quantities appear regularization independent. 

The partition function of the elastic material with one 
cut Zi in the saddle point approximation (^9|), will de- 
velop a multiplicative factor Zf upon inclusion of the 
quadratic fluctuations Zi ZjZi with 



Z/ = ^exp(-/3Ai?). 



(95) 



AS 



A deviation AE from the saddle point energy is decom- 
posed into three parts, with each part describing fluctu- 
ations of one of the mentioned three types 



n— 1 

inuous- 



1 + X 1 



(96) 



The first term in (|96| ) accounts for the decrease in the 
energy release due to the curving of the saddle point cut 
of length £c with the curvature 



/2A irnx 
k(x) = Z^c„W — sin-^, X € [OJc\- 

n—l ^ ^ 



(The first term in ( |96[ ) follows from (M) with / = Ic 
given by (^2|).) The second term in ( P6|) describes the 
asymmetric modes in the thermal fluctuations of the free 
surface of the saddle point crack about its equilibrium 
opening shape 



^asym^^^ = ^ (m„ + W„) siu 717?, 7? G [-TT, Tt) 



(97) 



'it) 



where a point at the cut boundary is parameterized 
by its distance t — £c{^ + cos7?)/2 from the cut end; 
7? £ [~TT, 0) parameterize the lower boundary displace- 
ments and 7? S [0, 7r) parameterize the displacements of 
the upper boundary points. The symmetric modes of 
the crack opening about its equilibrium opening shape 
are assumed to relax providing the minimum increase 
in the elastic energy for a given {7i„,7;„}. The latter 
guarantees that all additional modes with the continuous 
displacement at the cut boundary (the ones which give 
discontinuous — the last term in (^ describing the bulk 
phonons) decouple from {c„,?/„,7;„} and are the same 
as the ones for the uncracked material. (The arguments 
here are the same as those that were used in derivation of 
(^).) Since the curvature modes {cn} give the equilib- 
rium energy of the curvy cut, the response of the surface 
phonons to such a curving is already incorporated, so 
the quadratic fluctuations {cn} can be calculated inde- 
pendently from the quadratic fluctuations {?x„,7;„}. The 
latter means that there are no coupling between {c„} 
and {un, Vn} modes in (p6|), and the spectrum of {?/„, Vn} 
modes is the same as that for the straight cut of length 
4 (0). 

The last thing we have to settle before the calculation 
of Zf is the proper integration measure for the surface 
phonon modes {7i„,Wn}- (We argued in the conclusion 
of section III that the integration measure for the modes 
On is Cartesian — O^^i dcn-) Here we show that be- 
cause the functional measure in the displacement fields 
{u{x,y),v{x,y)) defined at each point of the material 
{x, y) is naturally Cartesian — D[u{x , y) / X]D[v{x , y) / X], 
the integration measure for the modes {7i„,7;„} must be 
of the form Yl'^^^{l/2n)dundvn/ ■ 

An arbitrary elastic displacement field for the mate- 
rial with a curvy cut is defined by specifying its bulk 
part (7ibuik(a;,y),7;buik(2^,y)) (point (x,y) can be any- 
where except at the cut boundary) and the cut part 
(w+t(*)>w™t(^):'fc^ut(i)>^'c"ut(*)) (the cut displacements 
are defined along the cut and are parameterized by the 
distance t = £c(l + cos7?)/2, 7? G [0, tt) from the cut end; 
the + and — superscripts are correspondingly the dis- 
placements at the upper and the lower boundary of the 
cut). It is helpful to visualize the introduction of the cut 
into the material as splitting in half each of the atoms of 
the material along the cut boundary. Then, the bulk part 
of the displacement field combines degrees of freedom of 
all atoms left untouched by splitting and the cut part 
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describes the displacements of the spht ones. Note that 
the sphtting increases the total number of the degrees of 
freedom. The original measure is naturally 

D [wbuik ix,y)/X]D [vhuik {x,y)/X]D [ucut {t)+ Ucut (ty /X^] 

D[Vcntit) + V,^t{t)-/X^] 

First we separate the symmetric and asymmetric parts 
in the crack opening displacement 



(99) 



Because the Jacobian of the transformation 
is constant 



the integration measure remains Cartesian: 



4, (100) 



^Kuik(a;, y)/X]D[v^M^, y) / XMu'^^ {t)v'y'^ (t) / X^] 

D[u''''y'^{t)v^'y"\t)/4X% 

Now we can combine the bulk and the symmetric cut 
part of the measure by introducing the continuous dis- 
placement fields (uc{x, y) , Vc{x, y)) everywhere, including 
the cut boundary. (In our atomic picture, the symmetric 
modes of the cut part of the displacement fields repre- 
sent the displacements of the split atoms as if they were 
whole, and so it is natural to combine these degrees of 
freedom with the bulk ones. Obtained as a result of such 
combination the continuous degrees of freedom are indis- 
tinguishable from the degrees of freedom of the uncracked 
material.) The integration measure becomes 

D[u,{x,y)/X]D[v,{x,y)/X]D[u^y''\t)v^y'^it)/4X% 

According to our decomposition, we specify the asym- 
metric cut opening and find the equilibrium displace- 
ment fields that minimize the increase in the clastic en- 
ergy. In other words, given {u'^^y"^ {t) , v^^^^"^ (t)) deter- 
mine (u™'" (u'^^y'" , v'^y"' ) , w™'" (u'^^y" , v'^y"' ) ) . The trans- 
formation 



u,ix, y) = ur\u^y''\ v^^"") + u,{x, y) 
Vcix, y) = vt^{u^^'^, v^^y^) + v,ix, y) 



(101) 



then completely decouple the surface phonon modes and 
the continuous modes that contribute to Ai^continuous in 
( p6| ) . The Jacobian of the transformation 

{u,{x,y),v,{x,y),u^^y^{t),v^^y^^{t))-^ 
{u,ix,y),Ux,y),u^y^^it),v^y^\t)) 

is unity (the transformation is just a functional shift) and 
so the measure remains unchanged 

D[u,{x,y)/X]DMx,y)/X]D[u-'y'^{t)v^y'^{t)/4X']. 



The Fourier transformation ([O^) is orthogonal, but the 
Fourier modes are not normalized: 



d-d sin n-d — —. 



(102) 



The latter means that at the final stage of the change 
of variables {u'^^y'^{t),v^^y'^{t)) {un,Vn} there appear 
the Jacobian n^i(2/'''')i ^nd so we end up with the in- 
tegration measure J^^j^(l/27r)du„du„/A^. 

From (|95|-p6|) with the proper integration measure over 
the surface phonon modes we find 



00 (.+00 f 



00 

n 

n=l 



Tr'^Xn'^ 



(103) 



+ °° 1 dUndVn , r, 
2 '^^P'l 



-00 



2^ A2 



-r—\Un+'"n) 



continuous 

3\„2\ 1/2 00 



n 



TT-^Xn^ 



n— 

11 zLtt/i 



-'continuous 



(104) 



where 



continuous 



E 



exp(-/3A£;continuous)- (105) 



A Bo, 



Because Ai?continuous corresponds to the degrees of free- 
dom of the uncracked material (with the same energy 
spectrum), ^continuous contributes to the partition func- 
tion Zq of the material without the crack, which accord- 
ing to (|9^) drops out from the calculation of the imagi- 
nary part of the free energy. 

All the products over n in these expressions diverge: 
we need a prescription for cutting off the modes at short 
wavelengths (an ultraviolet cutoff). 

First we'll consider the ^-function rcgularization. In 
this rcgularization prescription |^l| , the infinite product 
of the type D = is evaluated by introducing the 

function Dq{s) 



so that 



00 



Z? = exp(-i?'(0)). 



(106) 



(107) 
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It is assumed that the sum (106) is convergent in some 
region of the complex s plane and that it is possible to 
analyt ically continue Di^{s) from that region to s = 0. 
From (104) we find 



Zf = D1D2Z, 



continuous 



(108) 



where Di and D2 are obtained following (107) from the 
corresponding ('-functions: Dii^{s) and D2(^{s) 



Dicis) 



V 1 + X 



n— 1 ^ 



X 



Cr{s) in (109) is the standard Riemann (-function, holo- 
morphic everywhere in the complex s plane except at 
s = 1. Noting that ((0) = -1/2 and C'(0) = -(ln27r)/2 
we find from ([l07| - p^ 



WA j [l + x 



-'continuous- 



(110) 



From ( |92[ ) and (104) we find the imaginary part of the 
free energy in the (-function regularization 



J ] ImF'^""?''^ (Ill) 



where IniF'^™P''' is given by (|9S 

Second, we consider lattice regularization — which is 
more elaborate. We represent a curvy cut by TV -I- 1 = 
£c/A segments of equal length parameterized by the kink 
angles {ai}, i € [1,-^]- With our conventional param- 
eterization of the cut t — £c{l + cosi?)/2, tt, tt), the 
asymmetric modes of the crack opening displacements 
{u'^y"^ (t) , v^^"^ (t)} are linear piecewise approximation 
for given asymmetric displacements of the "split" kink 



atoms { 



asvm asvm 



^}, i e [1,^]- More precisely, if ti 
and ij+i parameterize the adjacent kinks, we assume 

asvm asym 

M^«y-(t) = u^^""'' + — {t - u) (112) 

asym asym 



V 



for t e [ti,tij^i\. From the integration measure ar- 
guments for the (-function regularization, it is clear 
that the integration measure in this case must be 



rii^i '^Q^i Dill du^^™ dv^^™ / 'iX^ . Now we have to write 
down the lattice regularization of the quadratic devia- 
tion A£' from the saddle point energy (96). From (^) 
and (^8|), the curving of the critical cut Ic will reduce the 
energy release by Ai?c 



32^ 

N 



N 



where 



From I 
by 



min(i, j) 



(113) 



(114) 



{N + lf N + 1 ' 
the surface phonon contribution to AE is given 



AEr. 



where from ( |9^ ) 



27r/i 



+ 00 



- / dd w^">""(t(z9))sinni? 



(115) 



(116) 



In principle, for a gi ven piecewise approximation of the 



asymmetric modes ( 112 ) determined by {u^^^"^ , vf^^"^} , 
i G [li^li o ne co uld calculate the Fourier amplitudes 
according to (116) and then plug the result into (|ll5| ) 



to obtain AEp in terms of {u'^^^'^ , v^^^"^} . We will use 
another approach. Using u^^y"'{0) = ■«^">"°(4) = (with 
the same equalities for v'^^'^''^ (t)) we integrate dll9) by 
parts to obtain 



dd 
d^ 



du^^y'°(i(i9)) cosn?? 



(117) 



Substituting (117) into (115) we find 



AEp = -jP^, I I d^id^2 

7^(1 +X) Jo Jo 



+ 



d^i 



ddi 



where 

00 

n=l 

Following H 



ECOS KX 1 , 
= - in 
k 2 

fc=i 



cosnt?i cos ni?2 



(118) 
(119) 



1 



2(1 — cos x) ■ 



(120) 



we find an analytical expression for the kernel (119) 
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K{'di,d2) = -^ln2- iln|cos?9i -cosi?2|. (121) 
Finally, introducing Aff from 



N 

E 



^asym^asym ^ ^asym^asym 



K{di,d2) (122) 



we obtain 



8/^ 



-(1 + x),^, 



E/ asvm asym , a 



asym asym , asym asym \ 71 



(123) 



To calculate we substitute (112) directly into the 
RHS of ( |122| ) and read off the corresponding coefficient, 
given by the following three equations: 

Mf, = h[h]) + h[i + + 1) - hit + 



-/2(*,J + 1) 
3 



4 (cosi?i_i — COSl?i)(cOSl?j_i — COSl?j) 
hih]) + h{i - 1, J - 1) - h{i - 1, j) - hih] - 1) 



/l(»,j) = 

(cos — COS i9j )^ 
4 

0, 



In 



sm — sm — 



if « 7^ j; 
otherwise, 

(124) 



where = 0; ^^Af+i = and 'di parameterizes the i-th 
kink (kinks are equally spaced in real space): 



= arccos 1 



2i 



N +1 



i G [l.,N]. 



(125) 



From ( 113 ) and (123) we find the quadratic deviation 
from the saddle point energy 



AE — AEc + AEp + Ai?continuous 



(126) 



N 



asym asym\^p 



asym asym 



continuous ■ 



Thus the multiplicative factor Zf to the partition func- 
tion of the elastic material with one cut in the lattice 
regularization is given by 



Zf^Wj danexpl -fiai^ ^ a.a-jM^j I (127) 
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exp< — p 



7r(l + x) 



asym,, asym _l asym asym 1 , y^p 



/ TT 



AT 



det -^continuous 



where ^continuous is given by ( |105| ). 

The determinant coming from the curvy cuts M^j can 
be calculated analytically. In section III we show that 
siuTrnx/^ are eigenvectors of the operator (^o|), the con- 
tinuous analog of M^j . One can explicitly check that for 
n G [1, A^], vectors m„ — {smTTni/{N + 1)} are in fact 
eigenvectors of M^j with eigenvalues 



1 



4(iV+l) 2{N + 1) 



and so 



N 



detM^^. = []A„= [N+l 



n=l 



-(N+l) 



(128) 



(129) 



(To obtain ( 129 ), we take the limit x ^ of 

27V+1 



sin2(7V-h l)a; = 2^^+! J| sin 
H, to get 



k=0 



2{N+1) 



(130) 



2N+1 



iV + 1 = 4^ Jl sin 



kn 



N 



k=l 



2(7V+1) 



4-n 



kn 



sm 



fc=i 



2{N + 1)' 
(131) 



With (131), the calculation in (129) becomes straightfor- 
ward.) Recallin g th at TV -1-1 = ^c/A, we can rewrite (127) 
making use of (129) 



(3aXf IT V"''^^ ( l^V'^ / 32/3/iA2 



TT \ I3a\ 

^^(1 + x) 

32/3/1A2 



7r2(l + x) 



£c/2A 



det^^Aff, Z, 



continuous- 



(132) 



Note that the first three factors on the RHS of (132) 
(coming from the curvy cut fluctuations) have the asymp- 
totic form, N ^ oo, 
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l(3aX 



1/2 



N'=^ exp{co + ciN} 

(133) 



with Co = 0, ci = \n{TT/(3aX)/2 and C2 = 1/2. 

We were unable to obtain an analytical expression for 
the surface phonon determinant det Mfj. For = 2...100 
kinks we calculate the determinant numerically and fit its 
logarithm with f{N) = po + piN + p2lnN (Figure 0), 



det AfP- = exp{po + Pi^V}. 



(134) 



We find Po = 0.09 ± 0.02, pi = 0.166 ± 0.002 and 
P2 — 0.24 ± 0.05. (We expect that the surface phonon 
fluctuations contribute to Zf similar to the curvy cut 
fluctuations (|T3^) — hence the form of the fitting curve 



for detAff .) 
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FIG. 7. We calculate numerically the logarithm of the sur- 
face phonon determinant, det M|^., for N = 2 ... 100 kinks and 
fit the result with /(TV) = po + P\N + p2 In 



Because we introduce new degrees of freedom with our 
"splitting atoms" model for the crack, we discuss the ef- 
fects of the corresponding new momenta. Before "split- 
ting", [tc/X — 1) atoms along the cut contribute 



^kinetic 
c/A-1 

n 



dpi^dpi. 



\2TTh 



exp< — /3 



^c/A-l 2 
PI 



E 



2ra 



(137) 



to the partition function of the uncracked material Zq. 
(We do not consider the contribution to the partition 
function from the bulk atoms — they contribute in a 
same way to Zq as they do to Zi , and thus drop out from 
the calculation of the imaginary part (p2|).) The config- 
uration space integration measure for a classical statis- 
tical system is dxdp/2'Kfi] because we integrated out the 
displacements with the weight 1/A to make them dimen- 
sionless, the momentum integrals have measure dpX/2Trh, 
( |137| ). The formation of the cut increases the number of 
the kinetic degrees of freedom by {£c/ A — 1) (the number 
of split atoms). The split atoms contribute 



kinetic 



n / / '^P'■.dp^y 



(—)' 

\2TTh J 



(138) 



to the partition function of the material with the cut. 
From ( |92| ) , ( |137| ) and ( |138| ) the kinetic energy of the elas- 
tic material modify the imaginary part of the free energy 
by a factor Z^: 



ImF 



lattice 



ZklmF 



lattice 



(139) 



From (132-134) 



Zf = exp{pi - Po} 



/32^VaA3^i/2^^^X-P2+i/2 



U'(l + X)y VA 
7r3(l + x)exp{-2pi}^'=/'^ 



32/3V"A3 



-'continuous 



(135) 



which following (|9^) gives the imaginary part of the free 
energy in the lattice regularization 



"P^^^"^"^^5(TT^ 
^3(1 -f;^) exp{-2pi} 
32/32^aA3 



with 



Zk = 



kinetic 
kinetic 



• 2 X 4/A-l 

mX ^ 



(140) 



One might notice that for both (^-function and lattice) 
regularizations the effect of the quadratic fiuctuations 
can be absorbed into the rcnormalization of the prefac- 
tor of the imaginary part of the free energy calculated in 
a simplified model (without the quadratic fiuctuations) 
(p3|), and the material surface tension a: the multiplica- 
tive factor to the imaginary part of the free energy has a 
generic form 



no 



X 



exp< n2 



X 



IwlF 



'Simple 



(141) 



with ImF'^'^^Pi^ from (|3|). 

Throughout the calculation we ignored the kinetic 
term in the energy of the elastic material: their behav- 
ior is pretty trivial, as momenta and positions decouple. 



where the first two terms renormalize the prefactor 
of lnii^''™pi'= and the other one can be absorbed into 
ImF^™P''' through the effective rcnormalization of the 
surface tension 
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2f3X 



n2 



(142) 



From (13£), (14C) it follows that in case of the lattice 
regularization, the inclusion of the kinetic energy of the 
elastic mater ial shifts the constants uq and n2, thus pre- 
serving (141): 



no no 



n2 



n2 + In 



m 



(143) 



The calculation of the kinetic terms in the C-function 
regularization is more complicated. We howev er ha ve no 
reason to believe that it will change the form (141). 



VI. 



THE ASYMPTOTIC BEHAVIOR OF THE 
INVERSE BULK MODULUS 



In our earlier work |13[ , we discussed how the ther- 
mal instability of elastic materials with respect to frac- 
ture under infinitesimal stretching load determines the 
asymptotic behavior of the high order elastic coefficients. 
Specifically, for the inverse bulk modulus K(P) in two di- 
mensions (material under compression) 



1 1 /dA\ 



c„P" 



(144) 



we found within linear elasticity and ignoring the 
quadratic fluctuations, 



Crt+l 



1/2 f Ax + I) 

64/3/ia2 



1/2 



as 71 ^ oo, 



(145) 



which indicates that the high-order terms Cn roughly 
grow as (n/2)! and so the perturbative expansion for the 
inverse bulk modulus is an asymptotic one. 

In this section we show that, except for the temper- 
ature dependent renor maliz ation of the surface tension 
a ar = a + 0(1//?), (145) remains true even if we in- 
clude the quadratic fluctuations around the s add le point 
(the critical crack); moreover, we argue that (145) is also 
unchanged by the nonlinear corrections to the linear elas- 
tic theory near the crack tips. 

We review how one can calculate the high order coeffi- 
cients of the inverse bulk modulus ||l^ . The free energy 
F{T) of the elastic material is presumably analytical in 
the complex T plane function for small T except for a 
branch cut T e [0, -l-oo) — the axis of stretching. (We 
show this explicitly in the calculation within linear elas- 
tic theory without the quadratic fiuctuations [|l3|.) It 
is assumed here that neither nonlinear effects near the 
crack tips nor the quadratic fluctuations change the an- 
alyticity domain of the free energy for reasonably small 
T (i.e. T <Y). One can then use Cauchy's theorem to 



express the free energy of the material under compression 
F{-P), (Figure I): 



n-p) = i 



2TTi J T + P' 

7 



dT. 



(146) 



A Im[T] 




FIG. 8. The free energy of the elastic material F{T) is 
analytical in the complex T plane except for a branch cut 
T £ [0, -|-oo). This allows a Cauchy representation for the 
free energy F{—P) of the material under compression. 



The contribution to (146) from the arc EFA goes to 
zero as the latter shrinks to a point. In this limit we 
have 



F{-P) = 



1 

2ni 



F{T + iO) - F{T - iO) 





1 
2^ 

B 



T 

F{T) 



P 



dT 



T + P 



dT 



BCD 



ImF(r) 
T + P 



dT 



1 

27ri 



F{T) 
T + P 



dT. (147) 
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As it was first est ablis hed for similar problems in field 
theory |2^ — ||2^, (147) determines the high-order terms 
in the expansion of the free energy F{—P) = /n-P" 



In 



(-1)" [ ImFjT) 



dT- 



(-1)" / F(T) , , 

BCD 



The second integral on the RHS of (148) produces a con- 
vergent series; and is hence unimportant to the asymp- 
totics: the radius of convergence by the ratio test is of 
the order the radius of the circle BCD (i.e. larger than P 
by construction) . The first integral generates the asymp- 
totic divergence of the inverse bulk modulus expansion: 



19 



(-1)" / Imf (T) 



dT as 



(149) 



Once a perturbative expansion for the free energy is 
known, one can calculate the power series expansion for 
the inverse bulk modulus using the thermodynamic rela- 
tion 



so that 



1 1 f dF(-P) 

K{P) ~ Tiy ap 



Cn+l (?l + 3)/„+3 



{n + 2)/, 



n+2 



(150) 



(151) 



Note that because the saddle point calculation becomes 
more and more accurate as T — s- 0, and because the in- 



tegrals in equation (149) are dominated by small T as 
n ^ oo, using the saddle-point form for the imaginary 
part of the free energy yields the correct n ^ oo asymp- 
totic behavior of the high-or der c oefficients /„ in the free 
energy. Following ( p3| ) and (141) the imaginary part of 
the free energy including the quadratic fluctuations is 
given by 



Imi^(r) 



1 



2rio 

32/iQ; 



nin2 \ ( 2PfiX'^ 



Til 



^r2(x + i)A 



V X 

A 
X 



1/2 



7rr2(x + l 

(152) 



where i s given by (142). Note that no, rii, n2 and 
ar in ( |l52| ) are regularization dependent coefficients, by 
our c alculations in the previous section. From (149) and 
(|l5l|) we find 



f Ti{x + l) \^'^ (n + 3)r(n/2 -h Til 2) 
V 32/3/ia2 



(n + 2)r(n/2-f ni -I-3/2)' 
(153) 



(In the limit n 
Using 



oo (153) is independent of B in (149) 
T{n/2 + ni + 2) 



r(n/2-Fni -H3/2) 



as n 



we conclude from ( |153| ) that 
64/3/iQ;2 



Cn+l 



as n - 



(154) 



(155) 



Equation (155) is a very powerful result: it shows that 
apart from the temperature dependent (regu lariz ation 
dependent) correction to the surface tension (142), the 
asymptotic ratio of the high order coefficient of the in- 
verse bulk modulus is unchanged by the inclusion of the 



quadratic fluctuations (at least for the regularizations we 
have tried). One would definitely expect the surface ten- 
sion to be regularization dependent: the energy to break 
an atomic bond explicitly depends on the ultraviolet 
(short scale) physics, which is excluded in the thermody- 
namic description of the system. This has analogies with 
calculations in field theory, where physical quantities cal- 
culated in different regularizations give the same answer 
when expressed in terms of the renormalized masses and 
charges of the particles |Q. Here only some physical 
quantities appear regularization independent. 

The analysis that leads to ( |155| ) is based on linear elas- 
tic theory that is known to predict unphysical singulari- 
ties near the crack tips. From ||l^, the stress tensor com- 
ponent (Jyy, for example, has a square root divergence 

(156) 



as one approaches the crack tip. One might expect 
that the proper non-linear description of the crack tips 
changes the asymptotic behavior of the high order elastic 
coefficients. We argue here that line ar an alysis gives how- 
ever the correct asymptotic ration (155): the /mear elas- 
tic behavior dominates the nonlinear asymptotics within 
our model. 

It is clear that the vital question is how the energy 
release of the saddle point (critical) crack is changed by 
nonlinear processes (microcracking, emission of disloca- 
tions, etcj in the vicinity of the crack tips as T ^ 0. Fol- 
lowing |2q ] we distinguish in the crack system two well- 
defined zones: the outer zone, consisting of exclusively 
linear elastic material, transmits the applied traction to 
the inner, crack tip zone where the nonlinear processes 
take place (Figure ^). Such separation introduces two 
length scales to the problem: rni and rcross- The first 
scale determines the size of the nonlinear process zone 
near the crack tips. It can be readily estimated from 
(156) by requiring the stresses at the boundary of the 



nonlinear zone to be of the order atomic ones ct, 



a 
Y' 



Y: 



(157) 



The second length scale is a crossover length rcross where 
the elastic fields near a crack tip deviate from the inner 
zone strain asymptotics to depend on the outer-zone 
boundary conditions (i.e. the length of the crack in our 
case). Normally, r™ss is only a few times smaller than 
the crack length [Q, — for the present calculation 
we assume rcross ~ ~ Ya/T'^, (p2|). 
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FIG. 9. In the crack system there are two well-defined 
zones: the outer zone, consisting of exclusively linear elastic 
material, and the inner, crack tip zone where the nonlinear 
processes take place. Such separation introduces two length 
scales: the first (rni) determines the size of the nonlinear zone, 
and the second (rcross) gives the scale where the elastic fields 
near the crack tip deviate from the ones predicted by SSY to 
comply with the outer zone boundary conditions. 

First, let's consider the energy in the nonlinear zone. 
The saddle point energy is aic and diverges as l/T^ as 
T — !■ 0, while the elastic energy in the nonlinear zone E^i 
is bounded by the linear value 



dr rajAr)/Y ^ a^/Y. 



(158) 



Since E-^i is fixed as T — > 0, it renormalizes no in (141) 
and hence does not affect the asymptotics (155). 

Second, we consider how the existence of the inner 
(nonlinear) zone changes the energy in the outer (lin- 
ear) zone. The elastic equations around the crack tip 
allow many solutions p^]; in each, the stresses cry have 



the form Cbr''fb{0), r„i < r < 



in polar coor- 



dinates (r, 0) centered at the crack tip, where b is an 
half-integer, the Cf, are constants, and the fb are known 
trigonometric functions. Linear fracture mechanics pre- 
dicts b = —1/2 to be the most singular solution (com- 
pare with (|l56| )) only because modes with b < —1/2 
would give rise to singular displacements at the crack 
tip. Incorporation of the nonlinear zone r < r^i however, 
removes this constraint. In other words, the nonlinear 
zone introduces new boundary conditions for linear elas- 
ticity solutions, allowing them to be more singular. The 
dominance of 6 = —1/2 solution is known as small scale 
yielding (SSY) approximation. Analyzing the mode III 
anti-plane shear fracture, Hui and Ruina argued |28| that 
SSY approximation becomes more and more accurate as 
e ~ ''ni/^'cross ^ 0. (They expect that the same result 
can be extended for mode I fracture.) Clearly, in our 
case e ^ as T ^ 0; thus the dominant contribution 
still comes from b = —1/2 solution. In fact, following 
p7| we expect 



E 



= {■■■, -7/2, -5/2, -3/2} 



E 



= { 1/2, 3/2, 5/2,---} 



(153) 



The inelastic stresses at the outer boundary of the non- 
linear zone r ~ r-a\ are of order y, thus from ( 159 ), for 
n < —1/2, Cn — 0(e^^/^) (recall that e = 7'ni/''cross ~ 
(T/Y)'^). These more singular terms in turn gener- 
ate corrections to C„ with n > 1/2 oi order 0(e). 
(One can see this from the fact that the dominant con- 
tribution from the more singular terms at r ~ ic i s 
C-3/2(^c/rni)"^/2 e.) The dependence of C„ in (|l59| ) 
on the polar angle 6 is implied. 

There is a formal analogy between the arguments pre- 
sented here for the stress fields in the crossover zone with 
the quantum mechanical problem of the bound states of 
the hydrogen atom. When we treat hydrogen nucleus as a 
point charge, for each orbital quantum number, the elec- 
tron wave function has two solutions near the origin (the 
position of the nucleus) : one is finite as r ^ and the 
other one is divergent . In a point charge problem 

one immediately discards the divergent solution because 
it can not be normalized and thus can not represent a 
bound state. However, in a finite-size nucleus model one 
notices that the electron wave function outside the nu- 
cleus is a mixture of the finite and the divergent solu- 
tions of the point charge problem. The normalization 
is resolved because inside the nucleus the electron wave 
function satisfies a different equation and becomes finite. 
The radius of the nucleus serves as a short-distance cutoff 
similar to r^i in the crack problem. 

The change in the contribution to the saddle point en- 
ergy from the outer zone as a result of the introduction 
of the nonlinear zone, 5i?outerj is given by 



dr r- 



,-{c„}-4{ci'"-} 

Y 



(160) 



The dominant contribution to (160) comes from the cross 
term between n = —1/2 and n = —3/2 corrections in 



(|159|) 



Y Y 



(161) 



the correction renormalizes therii coefficient in the imag- 
inary part of the free energy ( |l52| ) (regularization depen- 
dent in a first place), leaving the asymptotic ratio (155) 
intact. 

It is no surprise that the nonlinear effect s do not change 
the generic form of the imaginary part ( 152 ). The de- 
tailed nonlinear description of the crack tips is a specifi- 
cation of the ultraviolet (short scale) physics and thus is 
nothing but another choice of the regularization. From 
our experience with ^-function and the lattice regulariza- 
tions, we naturally expect that this nonlinear regulariza- 
tion preserves the form of the imaginary part (152). 
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Finally, let's consider the enhanced nucleation of sec- 
ondary cracks in the high-strain outer-zone region — a 
possible cause for breakdown of the "dilute gas" approx- 
imation. Inside the nonlinear zone of the saddle point 
crack, the critical crack length for a second crack is of 
the order ajY from ( ^2|) with T ^ Y ) and thus such 
microcracks can be easily created. In fact, the nucleation 
of these microcracks may well be the dominant mecha- 
nism of the main crack propagation. Micro-crack nucle- 
ation in the nonlinear zone will change the stress fields 
near the crack tips, but as we discuss above, has little 
impact on the saddle point energy (as the total energy in 
the nonlinear zone is finite) . We show now that such sec- 
ondary crack nucleation is exponentially confined to the 
nonlinear zone of the main crack. The probability Wiro) 
of the second crack nucleated somewhere at r > tq ~ r-^x 
(r = corresponds to a crack tip) is given by 



W^(ro) 



vdv 

^exp{-/3a^(r)} 



(162) 



where £{r) is a critical crack length at distance r from 
the tip of the critical crack. From (|2|) with T repl aced 
with the stress field near the crack tip given by (|l56| ), we 
find 



(163) 



(pa) with (163) gives 



+ CX) 

W{rQ) I ^ exp{-/3Q;r} = exp{-/3Q;ro} 
J {(iaX) 



(164) 



The exponential dependence of VF(ro) on t he bo undary 
of the nonlinear zone ro ~ r„; in equation ( 164 ) means 
that the nucleation of another crack (in addition to the 
saddle point one) is exponentially confined to the nonlin- 
ear zone, justifying the dilute gas approximation. 



VII. OTHER GEOMETRIES, STRESSES, AND 
FRACTURE MECHANISMS 

In this section we discuss generalizations of our model, 
more exactly its simplified version without the quadratic 
fiuctuations. We will do five things. In (A) we calcu- 
late the imaginary part of the free energy for arbitrary 
uniform loading and find the high-order nonlinear cor- 
rections to Young's modulus. We discuss the effects of 
dislocations and vacancy clusters (voids) in (B) and (C). 
Part (D) deals with three dimensional fracture through 
the nucleation of penny-shaped cracks: we calculate the 
imaginary part of the free energy and the asymptotic ra- 
tio of the successive coefficients of the inverse bulk modu- 
lus. Finally, in (E) we consider a non-perturbative effect: 



the vapor pressure of a solid gas of bits fractured from 
the crack surfaces, and show how it affects the saddle 
point calculation. 



A. Anisotropic uniform stress and the high order 
corrections to Young's modulus. 

We calculated the essential singularity of the free en- 
ergy at zero tension only for uniform isotropic loads 
at infinity. Within the approximation of ignoring the 
quadratic fluctuations, we can easily generalize to any 
uniform loading. In general, consider an infinite elastic 
material subject to a uniform asymptotic tension with 
CTyy — T, CTxx = eT (0 < e < 1) and a^y = 0. U sing the 
strain-stress analysis of ||l^ and following (^)-(|2l|), we 
find the energy E'reieasc, released from the creation of the 
straight cut of length £ tilted by angle from the x axis 



release 



64^ 



(l + e) + (l-e) cos 261 



(165) 



(The isotropic result (g) is restored for e = 1.) The new 
important feature that comes into play is that the crack 
rotation ceased to be a zero-restoring-force mode. Treat- 
ing the crack rotation to quadratic order in 9 from the 
saddle point value 6* = 0, we obtain the total energy of 
the crack E{M,e) similar to (|ll]8|), (|6|) 



E{M, 9) = a4 



(166) 



As before, Ai is the deviation of the crack length from the 
saddle point value £c, still given by (^2|). Following (p8|)- 
(p^), the imaginary part of the free energy for a dilute 
gas of straight cuts, excluding all quadratic fluctuations 
except for the uniform contraction-expansion (mode Ai) 
and the rotation (mode 9) of the critical droplet, is given 
by 



Imi^''™P''=(T,e) =± 



A 



1 



exp 



-32/3^a2 



1/2 



7rT2(x+l) 



(167) 



One immediately notices an intriguing fact: the e- 
dependence of the imaginary part is only in the prefactor, 
which, as we already know is regularization dependent 
anyway. In particular, the latter means that the inverse 
Young's modulus — the elastic coefficient corresponding 
to the transition with path e = — will have the same 
asy mpto tic behavior as that of the inverse bulk modu- 
lus (155): the asymptotic ratio of the high-order elastic 
coefficients of the inverse Young's modulus Y{P) 



1 



Y{P) A \dp)^ ° 



YiP----fy„P" + --- (168) 
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(P in (168) is a uniaxial compression) are given by 



Yr 



1/2 



as n — > oo. 



(169) 



B. Dislocations 

We have forbidden dislocation nucleation and plastic 
flow in our model. Dislocation emission is crucial for 
ductile fracture, but by restricting ourselves to a brittle 
fracture of defect-free materials we have escaped many 
complications. Dislocations are in principle important: 
the nucleation |^ barrier E^is for two edge dislocations 
in an isotropic linear-elastic material under a uniform 
tension T with equal and opposite Burger's vectors b is 



47r(l 



TTT In — 

2^ X 



En 



(170) 



where Eq is slT independent part that includes the dislo- 
cation core energy. The fact that E'dis grows like 1/ In T 
as T — > (much more slowly than the corresponding bar- 
rier for cracks) tells that in more realistic models disloca- 
tions and the resulting plastic flow cannot be ignored. 
While dislocations may not themselves lead to a catas- 
trophic instability in the theory (and thus to an imagi- 
nary part in the free energy?), they will strongly affect 
the dynamics of crack nucleation (e.g., crack nucleation 
on grain boundaries and dislocation tangles) |lj,p3|. 



C. Vacancy clusters 

We ignore void formation. It would seem natural to 
associate the negative pressure (tension) (— T) times the 
unit cell size with the chemical potential /i of a vacancy. 
At negative chemical potentials, the dominant fracture 
mechanism becomes the nucleation of vacancy clusters 
or voids (rather than Griffith- type microcracks), as noted 
by Golubovic and collaborators |Q. If we identify the 
chemical potential of a vacancy with — T, we find the 
total energy of creation a circular vacancy of radius R, 
^^vac(i?), to be 



(171) 



From (171) the radius of the critical vacancy is Rc — a /T 
and its energy is given by E^aciRc) = Tra^/T. A saddle 
point is a circular void because a circular void gains the 
most energy (~ area of the void) for a given perimeter 
length. In principle, the exact shape of the critical cluster 
is also affected by the elastic energy release. The latter, 
however, 



E, 
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8/x 



(172) 



is fixed as T — > 0, and thus the energy of the va canc y is 
dominated by £'vac(^c) for small T. (To obtain (172) we 
used the strain-stress analysis of [|l^ and expression ( |2l| ) 
for the energy release.) Using the framework developed 
for the crack nucleation, we find that in case of voids 
(again, ignoring the positive frequency quadratic fluctu- 
ations) the imaginary part of the free energy is given by 



Tm psiniplc (rp\ _ _)_ 



I (A 
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exp 
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(The special feature of the calculation ( |173| ) is that trans- 
lations are the only zero modes: the rotation of a circu- 
lar vacancy clus ter d oes not represent a n ew s tate o f the 
system.) From ( |173| ) we obtain following ( |l48|) and (151) 
the asymptotic ratio of the high-order coefficients of the 
inverse bulk modulus 



Cn+l 



(174) 



The divergence of the inverse bulk modulus is much 
stronger in this case: the high-order coefficients grow as 
c„ ~ n!, rather than as {n/2)\ (for the fracture through 
the crack nu cleat ion) . 

Whether (174) is a realistic result is an open ques- 
tion. Fracture through vacancy cluster nucleation is an 
unlikely mechanism for highly brittle materials: the iden- 
tification of with (— T) demands a mechanism for re- 
lieving elastic tension by the creation of vacancies. The 
only bulk mechanism for vacancy formation is dislocation 
climb, which must be excluded from consideration — the 
dislocations in highly brittle materials are immobile [p6| . 
Vacancy clusters might be important for the fracture of 
ductile (non-brittle) materials. However, the nucleation 
of vacancies must be considered in parallel with the nu- 
cleation of dislocations. Because at small T dislocations 



are nucleated much more easily ( 170 ) than vacancy clus- 
ters at low stresses, the dominant bulk mode of failure 
is much more likely to be crack nucleation at a disloca- 
tion tangle or grain boundary — as indeed is observed in 
practice. 



D. Three dimensional fracture 

Our theory can be extended to describe a three dimen- 
sional fracture transition as well. Studying elliptical cuts, 
Sih and Licbowitz |pl]] found that a penny-shaped cut in 
a three-dimensional elastic media subject to a uniform 
isotropic tension T relieves the most elastic energy for 
a given area of the cut. The energy to create a penny- 
shaped cut of radius R, £'pcnny(-R), is given by Q 



E. 



penny 



(i?) = 2TTaR^ - 



4(1 - a)R'^T^ 
3/i 



(175) 
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The zero modes contribute in this case a factor 27rV'/A^ 
— 2tt coming from the distinct rotations of the cut, and 
V/X^ coming from the translations of the cut. Here we 
find the imaginary part of the free energy to be 



1 



1/2 



exp 



3(1 -afr^ 



(176) 



and the asymptotic ratio of the high-order elastic coeffi- 
cients of the inverse bulk modulus 



/on n2 ^ 1/4 . .1/4 

Cn+1 ^ / 3(1 - g) \ / 



(177) 



From ( |l7q ) and (17£) it follows that the partition func- 
tion of the gas effectively increases the tension T by the 



vapor pressure Pyi 



r, where 
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/?A2 



vapor - -^cicp{-2l3aX). 



(180) 



In particular, the essential singularity of the free energy 
shifts from zero tension to minus the "vapor" pressure. 
This shift is clearly a nonperturbative effect. We were 
able to describe it only by allowing topologically differ- 
ent excitations in the system: a state of the elastic ma- 
terial with a bit completely detached from the crack sur- 
face may not be obtained by the continuous deformation 
of the crack surface (surface phonons) or the cut shape 
(curvy cuts). At zero external pressure, our material is in 
the gas (fractured) phase — not until Pvapor is the solid 
stable. 



E. Vapor pressure 

The approach we used to calculate the imaginary part 
of the free energy is a perturbative one. In a sense, noth- 
ing prohibits us from considering cubic, quartic, etc. de- 
viations from the saddle point energy. In fact, it is possi- 
ble to develop an analog of Feynman diagram technique 
(as in quantum electrodynamics |p9| or quantum field 
theory [^) and calculate the contribution to the imagi- 
nary part to any finite order. It is important to realize 
that even if we did this, the result would still be incom- 
plete: we would miss interesting and important physics 
coming from nonperturbative effects. Here we discuss one 
such nonperturbative effect, namely the "vapor" pressure 
of a solid gas of bits fractured from the crack surface. We 
find that including the vapor pressure, the essential sin- 
gularity shifts from T = to T = — Pvapor- Consider 
a dilute gas of straight cuts of arbitrary length with an 
elliptical opening (mode vi in (|7j)) and a solid gas of 
fractured bits from the crack surface. Following [Q, the 
partition function of the material with one cut Zi under 
a uniform isotropic tension T is 



Zi = ZqItt 
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1 



vi >Zgasie,vi), 



where Zgas(^,fi) is the partition function of the gas of 
fractured bits inside the crack of area tt£vi/2. It costs 
2aX to fracture one bit of size A x A from a crack step, 
so in an ideal gas approximation, the partition function 
of the gas is determined by 



Zgas(^,z'i) = £exp(-2/3aAiV) 



Tl = 



~2}^ 



expi -^jY exp(-2/3aA) 



VIII. SUMMARY 

In this paper we studied the stress-induced phase tran- 
sition of elastic materials under external stress: an elas- 
tic compression "phase" under positive pressure goes to 
a fractured "phase" under tension. We showed that in 
a properly formulated thermodynamic limit, the free en- 
ergy of an infinite elastic material with holes of prede- 
termined shapes is independent of the shape of the outer 
boundary as the latter goes to infinity. Under a stretch- 
ing load the free energy develops an imaginary part with 
an essential singularity at vanishing tension. To calcu- 
late the essential singularity of the free energy including 
quadratic fluctuations we determined the spectrum and 
normal modes of surface fluctuations of a straight cut, 
and proved that under the uniform isotropic tension a 
curvy cut releases the same elastic energy (to cubic or- 
der) as the straight one with the same end-points. The 
imaginary part of the free energy determines the asymp- 
totic behavior of the high-order nonlinear correction to 
the inverse bulk modulus |l3| . We find that although the 
prefactor and the renormalization of the surface tension 
are both regularization dependent (once we include the 
quadratic fluctuations) , the asymptotic ratio of the high- 
order successive coefficients of the inverse bulk modulus 
apparently is a regularization- independent result. 

Within our model, the asymptotic ratio is unchanged 
by the inclusion of the nonlinear effects near the crack 
tips. We generalize the simplified model (without the 
quadratic fluctuations) to anisotropic uniform strain and 
calculated the asymptotic behavior of the high order non- 
linear coefficients of the inverse Young's modulus. We 
computed the imaginary part of the free energy (and the 
corresponding divergence of the high-order coefficients of 
the inverse bulk modulus) for the fracture mechanism 
through void nucleation which dominates at small exter- 
nal pressures: we argue that it may not occur in brittle 
fracture and should be preempted by dislocation motion 
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in ductile fracture. We find that the simplified model ap- 
plied to three-dimensional fracture predicts a {n/4y. di- 
vergence of the nonlinear coefhcients of the inverse bulk 
modulus. 

Our results can be viewed as a straightforward exten- 
sion to the solid-gas sublimation point of Langer [ ^Jl9| 
and Fisher's |^ theory of the essential singularities at 
the liquid-gas transition. Indeed, if we allow for vapor 
pressure in our model, then our system will be in the 
gas phase at P = 0, as noted in section VII(E). The es- 
sential singularity we calculate shifts from P = to the 
vapor pressure. If we measure the nonlinear bulk modu- 
lus as an expansion about (say) atmospheric pressure, it 
should converge — but the radius of convergence would 
be bounded by the difference between the point of ex- 
pansion and the vapor pressure. 
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